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ABSTRACT 

We study conditions for formation of recollimation shocks in jets interacting with stellar winds 
in high-mass X-ray binaries. We show the existence of a critical jet power, dependent on 
the wind rate and velocity and the jet velocity, above which a recollimation shock is not 
formed. For jet powers below critical, we derive the location of the shock. We show that 
surface shocks may still exist above the critical power, but only occupy a small volume of 
the jet and do not significantly alter the jet opening angle. We test these prediction by 3-D 
numerical simulations, which confirm the existence and the value of the critical power. We 
apply our results to Cyg X-1 and Cyg X-3. 

Key words: acceleration of particles-ISM: jets and outflows-stars: individual: (Cyg X-1, 
Cyg X-3) - X-rays: binaries. 


1 INTRODUCTION 


Interaction of an extragalactic jet with a stellar wind of 
a massive star passing thr ough the jet was considered by 
iBednarek & Protherod l ll997h . As the shock-formation condition, 
they used the balance of the ram pressures of the forward motion 
of the jet matter and that of a radial stellar wind. This appears to be 
the proper approximate condition if the velocity of th e star within 
the jet is much less than the wind velocity. See also iKomissarovI 
( 1 19941) for the case of jet interaction with winds of low-mass stars. 

_ The same condition was then used bv IPubus, Cerutti & Henril 

( l2010h for formation of a recollimation/reconfinement shock in a 
jet originating from a compact object and interacting with the stel¬ 
lar wind of a massive do n or in a binary system. The condition for- 
mulated by iDubus et alj (1201 Oh was then applied to Cyg X-3 by 
IZdziarski et alj j2012bh . who concluded from it that a recollima¬ 
tion shock can be formed only if the initial jet half-opening angle, 
ag, is rather large, ^ 30°. 

However, the wind-jet collision in a binary system is clearly 
not head-on. Close to the binary plane, the wind velocity is parallel 
to the plane, and the jet is launched away from it. In the case of a jet 
perpendicular to the binary plane, see Fig. [T] the jet and wind ve¬ 
locities are close to perpendicular along the binary plane, and only 
a velocity component related to the jet opening angle is collinear 
with the wind. 


This effect was 
lytical estimates of 
after YH15 (see 


t aken into acco unt 
Yoon & Heinj j2015ll. 


here¬ 
sy and 


Perucho & Bosch-RamonI . 
[Perucho, Bosch-Ramon & Khangulvanll20ld for previous studies 
of formation of shocks in X-ray binary jets). As the condition of 


the ram pressure of the jet being in equilibrium with the wind ram 
pressure, YH15 took only the jet velocity component collinear with 
the wind direction in the vicinity of the binary plane, i.e., Vj sinoro, 
where vj is the jet velocity and ag is the initial opening angle of the 
jet. This is applicable only if the recollimation shock takes place at 
a distance from the jet origin of z <K a, where a is the separation 
between the binary components. They derived a formula for the 
location of the recollimation shock in this approximation. YH15 
also performed numerical simulations of the jet-wind interaction 
and found their formula agrees relatively well with the simulation 
results in the case of a moderate jet power, 10^^ erg s“', interacting 
with a strong wind with a high mass loss rate of = IO^^Mq 
yr“' and a high wind velocity of v„ = 2.5 x 10* cm s“'. However, 
the assumption that the shock takes place at z <K a is not applicable 
for a sufficiently strong jet. 

Here, we generalize the calculation of YH15 to the case of an 
arbitrary jet power. We concentrate on the c ase of systems in which 
the jet and orbital axes are aligned; see, e.g.. lSteiner & McClintockI 
( l2012h for a discussion of conditions for the spin-orbit alignment. 
We test the calculations by means of numerical simulations. We 
apply our results to jets of the high-mass X-ray binaries Cyg X-1 
and Cyg X-3. 


2 FORMATION OF RECOLLIMATION SHOCKS 
2.1 The setup of the prohlem 

Suppose the jet is expanding with an initial opening angle ag. The 
internal pressure, pj i oc (where y is the jet’s adiabatic index), 
of the jet drops more rapidly than the jet’s lateral ram pressure, 
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Figure 1. A sketch of a jet originating from a compact object (CO), such as 
a neutron star or a black hole, and a stellar wind. The jet direction is at an 
angle, with respect to the direction from the CO to the star. In the case 
considered here, the jet is along the z-direction, and thus f = 90°. The jet’s 
half-opening angle is a and the separation is a. The star, located at x, = a, 
emits an isotropic and uniform stellar wind. We also show a unit vector 
along the contact discontinuity, uq, and one perpendicular to it, uqj.. 


Pj,.r and the lateral expansion of the jet will become super¬ 

sonic. Whenever the internal pressure of the jet drops below the 
total ambient pressure, a shock will be driven laterally into the jet 
in order to bring the internal pressure into equilibrium with the bow 
shock. 

Such a shock will necessarily reduce the lateral expansion ve¬ 
locity of the jet (using the free energy in the jet’s expansion to in¬ 
crease the internal energy of the jet fluid), and may thus be called 
a “recollimation” shock. Below, we will distinguish “strong recol- 
limation shocks” from “weak” or “surface shocks”. A strong recol¬ 
limation shock will convert most or all of the lateral expansion en¬ 
ergy (a fraction ~ sin^ ag of the total kinetic energy of the jet) into 
particle heating and/or acceleration, with important observational 
consequences for, e.g., the high-energy emission from particles ac¬ 
celerated in this shock. A surface shock will reduce the opening 
angle by a small amount, Aa <K a, sufficient to raise the pressure to 
the ambient value, but dissipating only a small fraction of the jet’s 
lateral kinetic energy. 

To estimate the change in opening angle of the jet, we consider 
ram pressure balance of a supersonic spherical wind colliding with 
a supersonic conical jet. The collision will take place on the side 


of the jet exposed to the wind. At a given point along the jet, this 
can be approximated as two streams colliding at some angle. This 
collision will lead to formation of two shocks - one into the wind 
(which we will refer to as the bow shock) and one into the jet (the 
recollimation shock). 

We define the coordinate system with the x-y plane coinciding 
with the binary plane, and the z axis perpendicular to it, see Fig.[T] 
and, for the sake of simplicity, analyse the shock propeties in the 
x-z plane (where the bow shock and thus the recollimation shock 
will be strongest). We place the star centre at x, = +a and the jet 
origin at x = 0. Following YF115, we perform all calculations in the 
non-relativistic approximation. 

The contact discontinuity will be parallel to the post-shock 
velocities of the jet and the wind. In the frame of the contact dis- 
continuity, ram pressure ba lance takes the form of, e.g., eq. (11) of 
iNalewaiko & Siko^ j2009l) . 

PwK±=Pi^f±’ ( 1 ) 

where p„, pj are the densities of the wind and the jet, respectively, 
and v„j. and vj^ are the wind and jet velocity components perpen¬ 
dicular to the contact discontinuity, respectively. We can define a 
unit vector along the contact discontinuity and a vector perpendic¬ 
ular to it as 


Uo = (sino'', 0, cos a') = |-ro, 0, y 1 - x^j, 

Uo± = (-cosa',0, sina') = 0, xoj, (2) 

which defines xg and the inclination angle a' of the contact dis¬ 
continuity relative to the initial jet axis. Then v„i = v„ • Uqx and 
Vjj. = Vj • Uo±. We can then solve equation Q for xo(z). The trajec¬ 
tory assuming the post-shock jet fluid follows the contact disconti¬ 
nuity is then given by the solution of dx/dz = tan a' = Xgj yjl — x^. 
We note that we do not impose conservation of energy and momen¬ 
tum, assuming that energy can be dissipated and radiated and the 
wind flow can continue around the jet. Also, the above condition 
is not satisfied close to the jet origin because the perpendicular ve¬ 
locity components are not supersonic. We neglect this effect since 
we are concerned with regions where the velocity components are 
supersonic. 

The jet is initially freely expanding and conical with the ini¬ 
tial semi-opening angle, ag, related to its initial Mach number, Aij 
(YH15). The lateral ram pressure of the jet, i.e., along the x direc¬ 
tion, is 


Pj.x = Pj tan^ agvf = 


:'rz^Vj ’ 


(3) 


where we assumed Vj to be constant. In the second equality, we 
have expressed the pressure in terms of the kinetic power of the 
jet-l-counterjet in the rest mass motion, Pj = ;rpjVj^z‘ tan^ ag. Then, 
is independent of the jet opening angle. Note that we use the 
power of both jets here, which is of interest when considering the 
energy budget of the system, in particular comparing the jet and 
accretion powers. The same convention was used by YH15. 

On the other hand, the ram pressure of the wind in its motion 
component along the binary plane is 


Pw,x PwW.JC 


A^wVw / a^ 
4tra^ [a^ + ) 


(4) 


where M„ is the mass-loss rate in the wind and v„ is its velocity, 
assumed constant along the jet. The second factor on the right-hand 
side appears due to both the wind density and decreasing as 


© 2015 RAS, MNRAS 000,[T]{8] 





















Recollimation shocks in high-mass X-ray binaries 3 


+ z})- This factor is not present in eq. (8) of YH15 because 
they consider only shock formation close to the binary plane. How¬ 
ever, as we show below, recollimation shocks in jets of HMXBs can 
take place at distances comparable to the binary separation, or not 
happen at all, in which cases that term is important. 

We solve for the inclination angle of the contact discontinuity; 

j4Fj/(MwV„Vj) [tt^ + z^) - az 

tanar'=tanao—^-, (5) 

yj4Pj/ (MwV„Vj) (fl2 -I- z^) + z^ tan Uq 

with a deflection angle at the leading edge of the jet of Aor = ag-a'. 


2.2 Strong recollimation shocks 


The formation of a strong recollimation shock corresponds to the 
direction defined by equation (O being substantially different from 
the direction of the jet boundary facing the star, i.e., Ao- > a. 
The limiting case for the formation of a strong shock is thus 
a' = xo(z) = 0, which corresponds to the treatment given in section 
3.2 of YH15. 

Equating the pressures leads to a quadratic equation in (z/fl)^. 
The equation has real positive solutions for the jet kinetic power 
satisfying 


Pj < Per 


1 . 

TT^wV„Vj. 


( 6 ) 


Note that P^^ is independent of the jet opening angle. 

More generally, in the case of jets inclined with respect to the 
orbit, it is straightforward to show that the shock formation crite¬ 
rion is 


P < P 

^ j — ^ cr 


1 -t cos^\^ 
sin^ / 


(7) 


where ^ is the inclination angle between the jet and the orbital sep¬ 
aration vector, a. Note that the critical power is different for ap¬ 
proaching and receding half of an inclined jet, as the sign of cos f 
changes from positive to negative as a jet tilts further than 90° away 
from the companion, because the wind ram pressure decreases and 
lower powers are sufficient to suppress the formation of a recolli¬ 
mation shock. 

If the condition in equation ([6]l is satisfied, the jet lateral ram 
pressure, decreasing as z“^, becomes equal to the corresponding 
wind ram pressure at some point zi along the jet, and a strong rec¬ 
ollimation shock forms. Thus, the lower solution corresponds to the 
condition of shock formation. The shock location is then given by 


(^j^=(5-l- V(52 - 2(5, 

l/(2d), d » 1, 


d = 


M„VwVj 

~^P~ 


> 2 , 


( 8 ) 

(9) 


where the approximate equation 12) corresponds to the solution 
given by YH15 in their eq. (9). At the critical point of Pj = Pa, 
there is one solution at zi = a. Thus, under the adopted approxima¬ 
tions, a strong recollimation shock never forms at z > a. If Pj > Pa, 
the jet lateral ram pressure will be always greater than that of the 
wind. 

Once the shock takes place, the jet is no longer conical and 
freely expanding [such that equation (O no longer applies], but 
instead its thickness follows from pressure equilibrium between 
the jet and the wind bow shock, see section 4.1 of YH15. How¬ 
ever, equations US still remain valid. Thus, the jet becomes 
again freely expanding above Z 2 given by the upper solution of the 
quadratic equation obtained by equating these two equations. 


= 5 - 1 -H V(52 - 2d, d>2. 

a / 


( 10 ) 


2.3 Surface shocks 


Even if the strong shock criterion is not satisfied, a shock may still 
form if the internal pressure of the jet drops below the bow shock 
pressure. In that case, the maximum deflection angle, Acr, from 
equation l|5]l is small compared to a. While such a surface shock 
will necessarily be highly oblique and weak compared to a recol¬ 
limation shock, it is important to note that it may still be a strong 
shock in the sense of having a large Mach number and changing 
the entropy of the gas. 

In the limit of Aa <K Uq, the location of the maximum relative 
dynamic pressure, where Aa is maximized, is 


Zmaz ~ - (1 -I- sinao) a, (11) 

cos tto 

and, to the lowest order, the maximum deflection angle becomes 



valid in the regime Pj » Mv„Vj (and assuming the small angle 
approximation for ag). 

The pressure behind such a surface shock is still given by the 
bow shock pressure, Pshock ~ jWv^/(16to^) (the same as in the case 
of a strong recollimation shock). The volume occupied by the shock 
is roughly 


Itto? 


2n 


^ -agAa ^ yttofl 



(13) 


This situation is illustrated in Fig.|2l which shows the structure of 
the bow shock and the moderate oblique shock driven into the jet. 

Since we would expect the non-thermal emission from such 
a shock to be proportional to the volume of the shock and some 
function of the pressure, the total non-thermal emission from the 

shock should be proportional to -yjMv^Vj/Pj, i.e., it decreases with 
the square root of the jet power for super-critical jets that do not 
form strong recollimation shocks and for otherwise identical jet and 
wind parameters. 


3 COMPARISON WITH SIMULATIONS 
3.1 Methods 

In order to test the predictions for the formation of recollima¬ 
tion shocks, we carried out 3-D hydrodynamic simulations us- 
ing ELASH 3.3 g rid-based adaptive-mesh hydrodynamics code 
dErvxell et alJ2000h for testing the conditions in which a recollima¬ 
tion shock forms, using the same method as that of YH15. In YH15, 
the jet was modelled by a cylindrical nozzle with inflow-boundary 
conditions at the surface, injecting a bipolar outflow with a pre¬ 
scribed energy, mass, and momentum flux. This cylindrical nozzle 
was required to set a slow lateral outflow from the side walls of the 
cylinder due to the numerical stability. The injected materials from 
the lateral outflow interact with the stellar wind, losing the momen¬ 
tum and falling back to the jet, which yields artificial mass-loading 
on the jet. In the jet model with Pj < 10'^’ ergs^', the amount of 
the mass-loading and its elfect on the jet bending are negligible. 
However, as the jet kinetic power increases, the mass-loading be¬ 
comes noticeable and this effect should be taken into consideration. 
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Figure 2. Left panel: Visualization of the P3e37 simulation. Red contour surface: section of the jet boundary; blue contour surface: section of the bow shock; 
yellow contour surface: surface of the star; bottom surface: slice through the temperature of the stellar wind and bow shock; vectors: stream lines sampling the 
jet and the stellar wind/bow shock. The axes are in units of 10^^ cm. Right panel: slice through P3e37 simulation aX z = 1.17 x 10^^ cm showing the structure 
of the bow shock and the formation of a surface shock on the leading edge of the jet; top part: five pressure contours linearly spaced between 13 and 80 per 
cent of the theoretical stagnation point pressure of the bow shock; overlaid in yellow is the approximate location of the contact discontinuity; bottom part: 2-D 
velocity streamlines in the x-y plane showing the deflection of the wind in the bow shock and the acceleration of jet material away from the stagnation point 
behind the surface shock and in the expansion wave in the down-wind region of the jet. 
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Figure 3. The density contour maps in x-z plane for the Ple37 (left), P3e37 (middle) and Ple38 (right) models. The Ple37 model is from YH15, in which the 
considered length along the jet is an order of magnitude smaller than that of the P3e37 and Ple38 models. The vertical line up from the black hole position is 
shown by yellow dashes. The thick yellow line represents the identified jet centre from simulations and the thick red line represents the analytic approximation 
to the jet trajectory. 
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Figure 4. The jet in the y-direction as a function of z. The black solid line represents the jet identified from the simulations. The straight line between the black 
hole and the end of the jet is shown by the red dashes. In the Ple37 model, the blue dashes are straight lines that follow the initial opening angle, ag ^ 6.5°, 
and the horizontal dotted line represents the location of the recollimation shock. 


To eliminate the side effects of the artificial outflow, we implement 
the jet nozzle within the lower boundary instead of an interior grid. 

Our numerical co n figurat ion is similar to t hat i n 
IPerucho & Bosch-RamonI ( 1200 Sll and IPerucho et all dlOlOh . 
which showed jet-wind interaction in massive X-ray binaries by 
using 2-D and 3-D simulations, respectively. They concluded that 
jet is likely bent and disrupted for jet kinetic power < 10^^ erg s“' 
by the ram pressure of stellar wind from nearby O-type star, which 
is also consistent with YH15. 

The parameters of the jet and stellar wind are the same as in 
YH15, except for the jet kinetic power. Namely, we use I'j = 3 x 
lO^cms^*, v„ = 2.5 X 10*cms“‘, = 10“^ Mq yr“*, a = 3 x 

10*^ cm, and the Mach number of Afj = 30, which corresponds 
to one of the two recollimation-shock cases studied by YH15. For 
the assumed vj, v„, and the critical jet power, equation is 
Pa - 3.0 X 10^’ erg s“'. 

In addition to the recollimation case considered by YH15, 
with a kinetic power of Pj = lO^^ergs^' (Ple37), we have car¬ 
ried out two sets of runs with Pj = 3 x 10^^ and 10^^* erg s“*, which 
we denote as P3e37 and Ple38, respectively. In our coordinate sys¬ 
tem, the jet origin and the star are located at Xco = 0 and x, = a, 
respectively. A volume rendering and contour plot of the P3e37 
simulations is shown in Fig.|^for illustration. 


3.2 Jet bending 

In the presence of a stellar wind, the jet is bent by a transverse pres¬ 
sure gradient that the wind drives around the jet (see Fig. H. The 
bending has to be taken into account for identification of a recol¬ 
limation shock. Therefore, we identify the jet centre in order to be 
able to estimate the jet thickness in the y-direction, see YH15. Fig. 
[3] shows how the jet is deviating from its initial direction along the 
z axis. The yellow thick lines show the jet centre identified in the 
simulation data. We ensure that the simulations have evolved for a 
sufficiently long time for the initial bow shock of the jet to prop¬ 
agate well beyond the region of interest, allowing for interaction 
directly with the stellar wind. 

For the analysis of jet-bending, we assume the jet remains con¬ 
ical with a fixed jet opening angle, ag. This is in contrast to the 
analytic jet model in YH15, in which the jet is confined by pres¬ 


sure equilibrium between the jet and the wind (red solid line of the 
Ple37 model in Fig. O, and thus the jet is not conical. The for¬ 
mer is appropriate for powerful jets and prior to the formation of a 
strong recollimation shock, the latter is appropriate for weaker and 
strongly bent jets beyond the location of the recollimation shock. 
As we will show in Section [33] below, the more powerful jets in 
our case study (the P3e37 and Ple38 models) are not substantially 
recollimated, and have approximately constant opening angle, ag, 
so we adopt a conical shape for the analytic jet model. 

The derivation is similar to that of YH15. Assuming that the 
longitudinal jet velocity is constant and AIj » 1, the longitudinal 
jet momentum per unit jet length can be calculated by 

Oj = fdA^PjVj = (14) 

J vf 

where dAj^ is the area perpendicular to the initial jet direction. The 
transverse wind momentum per unit jet length accumulated by the 
jet as a function of z is 

i r dz'v^^(z')Pw(z')^(z'), (15) 

''j Jo 

where h(z) = 2z tan ag is the jet width. Equation {H can be readily 
integrated to 




v„,oM„tanao 

47rvj + zj 


(16) 


In the first-order approximation, the jet bending angle, (/r, can 
be derived by the ratio between the accumulated transverse mo¬ 
mentum and the longitudinal momentum as a function of z, 


l/r(z) 


Ad»» 


Vw.oVjMwtanor 7 } 
47rPj cP + 7 } 


(17) 


At ^ a^, the bending angle reaches its asymptotic value. In the 
small-angle approximation, we can derive the analytic jet trajectory 
through d(x - xco)/dz = tan {j/. Analytic trajectories derived for the 
simulations are shown by the thick red lines in Fig. They show 
good agreement with the simulation results (thick yellow lines), 
indicating that our simple analytic model can describe the bend¬ 
ing jet. The bending angle is similar between the P3e37 and Ple38 
models, because although the Ple38 model has a jet kinetic power 
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about 3.3 times larger than that of the P3e37 one, the opening an¬ 
gle in PIe38 is also 2.4 times larger than that of P3e37 (see Section 
[33}. 


3.3 Recollimation 

YH15 studied the appearance of a recollimation shock at Pj = 
10^’ergs^*, in which case the shock takes place at z 7 x 
10" cm <K a (horizontal dotted line in the left panel of Fig.|4). 

In YH15, the jet was found to initially expand freely, its lat¬ 
eral expansion proceeding until the lateral ram pressure of the jet 
dropped below the wind ram pressure, where a recollimation shock 
occurred (independent of the initial jet opening angle and the jet 
Mach number). Beyond the shock, the jet evolved adiabatically in 
pressure equilibrium between the jet and the wind bow shock. 

However, with increasing jet power, the shock takes place at 
increasingly larger z or does not occur at all, as found analytically 
in Section[2] To check for the presence of a recollimation shock, we 
measured the jet width facing the star, i.e., along the y-direction. We 
measured the width along the identified jet centre (yellow thick line 
in Fig.|3}. The black solid lines in Fig. [4] show that the measured 
jet width increases almost linearly, and no recollimation shock is 
apparent. The opening angles are Q'(P3e37) =: 3.4° and Q'(Ple38) - 
8 °. 

In order to determine whether the contact discontinuity is mea¬ 
surably misaligned with the initial jet opening angle, we compare 
the jet trajectories with straight lines connecting the black hole with 
the edge of the jet at the upper limit of the analysis box (red dashes 
in Fig.|4}. Then we calculate the (squared) fractional deviation as a 
function of z, 

o'^(z) = [yjb(z) - (18) 

where yj and yjo are the positions of jet edge for the simulations 
and the straight line, respectively. Fig.j^shows that for Ple38, the 
identified jet boundary is in good agreement with the straight line, 
indicating a conical shape of jet and the absence of a recollimation 
shock. 

The jet in P3e37 also has an overall almost conical shape, but 
there is substantial negative deviation from this shape at z ~ 3 x 
10" cm. This shows that a surface shock does take place around 
z ~ fl, confirming the analytic prediction of Section|2]for the shock 
location around the separation at Pj > Per- The jet at z ~ a has then 
an opening angle a > 3.4°. 

These results are in an good agreement with the predictions 
of equation ll} for the wind properties that we have used in this 
work. According to equation H}, the critical jet power, above 
which there should be no strong recollimation shocks, equals Pj 
3.0 X 10^^ ergs“‘. From our simulations, there is a strong shock at 
Pj = 10^’ ergs^', a mild shock at Pj = 3 x 10^^^ erg s“', and and no 
measureable change in opening angle for Pj = 10^^* erg s“'. 

While the analysis presented in this paper focuses on the shock 
formed by the interaction of the leading edge of the jet with the 
stellar wind bow shock, the simulations allow a description of the 
down-wind half of the jet dynamics as well. As can be seen in 
Fig .|2 the thermal pressure of the bow shock drops to low values in 
the shadow of the jet, leading to the formation of an expansion fan 
that accelerates the jet and wind fluid into the evacuated region be¬ 
hind the jet. In some instances, this very low-density material may 
form a shock upon convergence on the x-z plane behind the jet, 
which is the reason for the appearance of denser gas in the wake of 
the jet in Fig.|3 



Figure 5. The squai'ed fractional deviation between the jet edge measured 
in simulations in the y-direction and the straight line connecting the black 
hole and the jet end, given by cr^ defined by equation d). The black and 
red colours correspond to the P3e37 and Ple38 models, respectively. The 
horizontal dashed line indicates no deviation. 


4 JETS IN CYG X-1 AND CYG X-3 


We can check the implications of the above results for two high- 
mass X-ray binaries wi th jets. Cyg X-1 contains a black hol e with 
the mass of =: 16 Mq dOrosz et al.l 201 1|: Zi6tko wskill2014l) . The 
mass of the donor is disputed. Iziolkowskil dTOldll estimates it as 
27 Mq. Then, the separation equals 3.2 x 10" cm. Cyg X-1 has 
a steady compact jet in the hard state dStirling et al.l2001 ). 

The current best estimates of the relevant parameters of Cyg 
X-1 are the mass-loss rate in the hard state of =: 2.6 x 10 “'’Mq 
yr“‘ doles et al.l 20031. it s asym ptotic velocity of v„ =: 1.6 x 10® 
cm s * doles & Bolton! 198^ and the r e latively uncertain jet 
velocity of Vj ^ 0.6c dStirling et ^ l200ll: lOleissner et Zll2004 
iMalzac. Belmont & FabiarJ ~ 2009ll . These parameters corre¬ 


spond to the critical power of — 3 x 10®^ erg s“*. This 
is within the estimated Pj == (1-7) x 10®^ erg s“' (for the 
jet-fcou nterjet) from th e optical nebula presumably powered by 
the je t dSell et al.l 1201 5l : see also lOallo et ^ l2005l : [Russell et al.l 
l2007h . Thus, a recollimation shock can occur in the hard-state 
jet of Cyg X-1. The location of the shock is given by equation 
At Pj = 10®^ erg s“', the shock takes place at z = 0.32a. 
The jet power of Cyg X-1 has also been estimated in severa l 
theoretical models of Zdziarski^ Lubihski & Sikot^ 1 2012ak 


[Malyshev. Zdziarski & Chemvakovci d2013[) . [Zdriarski et al.[ 

d2014[) and YH15. For almost all of the model power was <K P^ 
(the only exception was one model in [Malyshev et al.l|2013L which 
those authors themselves considered unlikely). 

Although it is likely that Cyg X-3 also contains a black hole, 
the curre nt solution for its mass also allows the presenc e of a neu¬ 
tron star dZdziarski. Mikola iewska_& Belczytis y[2013h . We adopt 
7 X lQ“^Mfri vr“' (see [Zdziarski et alj[2013l and references 
therein), v„ ^ 1.6x 10® cm s“‘ dvan K erkwii k et al .Il996[) . Vj =: 0.5c 
dZdziarski et aD [2012bl:lDubus et al.[[201(lh and a = 2.4 x lO" cm 
(assuming the total mass of 13Mi^. lzdziarski et alj[2013h . We then 
obtain — 6.6 x 10®’ erg s“'. This is sim i lar to, o r somewhat be¬ 
low, the power estimate of [Zdziarski et alj 1 201_2b}) during periods 
of strong y-ray emission in the soft state dFermi LAT Collaboratiord 
[2009l : lThvani et al1[2009[) . based on the y-ray luminosity. 

Thus, a strong recollimation shock may occur in the hard-state 
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jet of Cyg X-1 and the soft state of Cyg X-3 during periods of strong 
y-ray emission. Since ~ Pa in both cases, the shock location is 
close to the maximum po ssible position of z ~ a. This, in fact, 
agrees w ith the analyses of iDubus et al.l l l2010tl and IZdziarski et all 
( l2012bll for Cyg X-3, based on modelling of the y-ray emission. 

In Cyg X-1, the formation of a recollimation shock at z 
a may thus explain an apparent d eparture of the structure o f 
the jet from the s tandard model o f iBlandford & Kbnigll lll979l) . 
see (lleinj l l2006l) . IZdziarskH ( 12012h and Zdziarski et al.l l l2012ah . 
Nam ely, the observed re solved fraction of the 8.4-GHz jet emis¬ 
sion dStirling et alj|200lh implies the location of the onset of the 
optica lly-thin synchrotron radio emission at z ~ 10*"' cm dHeinzI 
l2006h . which appears not compatible with the observed orbital 
modulation of th e 15-GHz emission having a depth of =: 30 per cent 
(Izd ziarskj|2()I^ . This modulation results from free-free absorption 
in the wind, and its large depth requires a substantial difference be¬ 
tween the paths through the wind between the superior and inferior 
conjunctions, which, in turn, im plies that a large part of the radio 
emission originates from z ~ a ( IZdziarskill2012h . A recollimation 
shock at z i a can lead to strong dissipation and radio emission 
(spatially unresolved) on top of the remainder of the jet emission, 
explaining the above discrepancy. 

On the other hand, we note that Cyg X-1 in its hard state is 
much less radio loud than Cyg X-3 in its hard state. For the ra¬ 
dio loudness defined as the ratio of vLy at 15 GHz to the bolo- 
metric luminosity (dominated by X-rays), both in units of the Ed- 
dington luminosity, this strik i ng dif ference is illustrated in fig. 7 of 
IZdziarski. Segreto & Poolevl ( l2016h . This ratio is ~ 3 x 10“* and 
~ 10“® in Cyg X-1 and Cyg X-3, respectively. This difference by 
a factor of ~30 could be explained by the shock in the former and 
the latter being weak and strong, respectively. 

Finally, we stress that bri ght radio jets a re common also in 
low-mass X-ray binaries, e.g., iFendeJ d2006l) . where no stellar 
wind-jet interaction occurs. Thus, the formation of a recollimation 
shock, either strong or weak, is clearly not a condition for a jet ra¬ 
dio emission; however, it is likely to modify the radio emissivity 
profile with respect to that corresponding to the absence of such a 
shock. 


5 CONCLUSIONS 

We have derived an analytical criterion for the existence of recolli¬ 
mation shocks in jets in high-mass X-ray binaries. We have tested 
it with 3-D numerical simulations and have found an good agree¬ 
ment with the analytic predictions. According to this criterion, a 
recollimation shock in jets of high-mass X-ray binaries can exist 
only below the critical jet power proportional to the mass-loss rate, 
given by equation (|6j. Below that critical jet power, the position of 
the shock is given by equation (ID. Both the condition for shock 
formation and the shock location are independent of the initial jet 
opening angle. Above the critical power, a highly oblique surface 
shock may still form to raise the jet’s internal pressure to the bow 
shock pressure, but such a shock will only change the opening an¬ 
gle by a small amount and will only occupy a small fraction of the 
jet’s volume. 

We have applied our results to the high-mass binaries Cyg X- 
1 and Cyg X-3. The presence of a recollimation shock is likely in 
both of those objects. In the latter system, we can identify the po¬ 
sition of the shock with the region of y-ray emission. The large 
difference in the radio loudness between those objects can be ex¬ 


plained by the presence of a weak and strong shock in Cyg X-1 and 
Cyg X-3, respectively. 
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